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Abstract 

A relativistic wave equation for spin 1/2 particles in the Melvin space-time, a space-time where 
the metric is determined by a magnetic field, is obtained. The energy levels for these particles 
are obtained as functions of the magnetic field and compared with the ones calculated with the 
Dirac equation in the flat Minkowski space-time. The numeric values for some magnetic fields of 
interest are shown. With these results, the effects of very intense magnetic fields in the energy 
levels, as intense as the ones expected to be produced in magnetars or in ultra-relativistic heavy-ion 
collisions, are investigated. 


I. INTRODUCTION 


In the past few years systems with extreme magnetic helds {B > lO^^G) have been 
proposed to exist. In the magnetar analysis, for example, magnetic helds of the order of 


at the surface 


fl. s 


and at the center are expected to exist. In ultra-relativistic 


heavy ion collisions, at ■J^= 200 GeV (at RHIC), the magnetic held is expected to reach 


values as high as lO^^G [sl-j^ and at LHC, at \fs= 7 TeV, 10^°G. 

When studying particles inside this systems (eletrons inside stars, and so on), the usual 
way is to solve the Dirac equation 6|, [7|, for example, and hnd the energy levels, and in 
general, the Minkowski space-time is considered. A question that is quite intriguing, is if 
when the magnetic held reach these values, some ehect of the structure of the space-time 
may be observed. This is the purpose of this work. 


Including the magnetic held in 


exist, as the Melvin metric 


ic jsl- 12 


he metric is not a trivial question but some solutions 


12| . where a magnetic universe, with a magnetic held in 
the 2 ; direction, is considered, or the Gutsunaev solution jl3|, [l4|, for a magnetic dipole. 
Of course it is always possible to study an arbitrary shape of the magnetic held and solve 
(at least numerically) the resulting Einstein equations, but with the objective of hnding 
analytical results, and with this procedure, exploring this ehect in a hrst approximation, we 
will consider the Melvin metric and hnd the wave equations for Dirac particles subjected to 
a magnetic held inside this metric. 

So, this paper will show the following contents: In section II, a brief review of the formu¬ 
lation of the Dirac equation in curved spaces will be made. In section III the wave equation 
in the Melvin metric will be worked out. The results and conclusons of this work will be 
shown in section IV. 


II. DIRAC EQUATION IN CURVED SPACES 

In this section a brief review about the wave equation for spin 1/2 particles in curved 
spaces will be made. The basic formulation and the equations that will be needed in the 
next sections will be shown. 

A fundamental characteristic of the Dirac equation is its invariance under Lorentz trans¬ 
formations, so, when studying Dirac particles in curved spaces, it is interesting to preserve 









this aspect. One way to do that is by using the tetrads that may be dehned in order to 
satisfy the equation 


= e^%e^%Via)ib}, ( 1 ) 

where ?7(a)(6) is the Minkowski tensor, that represents a flat space-time, and the metric 
tensor related to a space-time that possesses an arbitrary geometry l^-[l7 . 

Observing eq. CD. we may note that the tetrads may be used in order to project vectors 
from the curved space-time in the flat space-time with the expression 


K = 


(a), 


( 2 ) 


that relates the form of a vector in different space-time geometries. 

As it was said before, the Dirac equation when written in the Minkowski geometry pos- 
seses Lorentz symmetry, and we want that the equivalent equation, written in a curved 
space-time possesses the same characteristic. In fact, if we study the behavior of the ele¬ 
ments of the Dirac equation under transformations that preserve the desired symmetries, we 
may understand which are the quantities that need to be added or modihed in the equation. 

We must note that a spinor transforms according to 


(3) 

where p (A) = 1 -|- and S(a)(6) is the spinorial representation of the generators 

of the Lorentz transformation, written in terms of the matrices, '^{a){b) = [ 7 (a )5 7(fe)] 

|i^ . The idea is to construct a covariant derivative V(a)'0 that is locally Lorentz invariant, 
that means that we need to impose the transformation condition 

^ pA (4) 

The usual way to obtain the form of the covariant derivative operator of the spinor is by 
supposing the combination 

= + (5) 

and considering the form of the operator transformation 

-)■ -h df,pp~^. (6) 

If we choose the combination of terms 


m 




or in an equivalent form 


O = -7T / N 

— 2*'^ (a)M(&)^ ’ 


( 8 ) 


with the term T (a)fi{b) dehned as 


r {a)fj.{b) — e(a)!7 + r ) ) 


(9) 


where are the Christoffel symbols, eq. (E]) is satished. Consequentlly eq. ([5]) satishes 
the transformation condition demanded in (jl]) and then we get the final form of the covariant 
derivative operator 

V(e)V' = e(^) ( 10 ) 

We must remark that there are different definitions of covariant derivatives of tensors in the 
literature when studying the covariant Dirac equation, as for example in 
Considering the Dirac equation in a flat space-time 


19|- 


2 l| 


5(a) "0 — ■m.'*/’ = 0, 


( 11 ) 


and replacing the conventional derivative operator by the one obtained in ffTOj) we obtain 
the desired wave equation for spin 1/2 particles in a curved space-time 

(5^ + 0.^,)^) -m'lp = 0 . ( 12 ) 

If the particle is submited to an external electromagnetic field, we may introduce this effect 
by a minimal coupling 


(^ i * + 'ill-mi’^ 0 , 


(13) 


and the scalar action that leads to this equation is given by 


S= d^x^J—gip 


+ leA^ ) - m 


ij. (14) 


It is usual to define the term 7 ^ = as a Dirac matrix in a given curved space-time 

and it is easy to verify that it satisfies the Clifford algebra 


YY + YY = 


( 15 ) 







III. WAVE EQUATION IN THE MELVIN METRIC 


If now we want to obtain the wave eqnations shown in the last section, in a space¬ 
time that has its strnctnre determined by a magnetic held, an nsefnl idea wolud be to 


consider the Melvin metric 


a- 


12| . The Melvin metric is a solntion of the Einstein-Maxwell 


^nations of the general relativity that represents a cylindrical magnetic nniverse. In his work 
8|, {q], Melvin considered a static magnetic held where its lines lie in cylindrical snrfaces 
perpendicular to the radial direction, with intensity ~ Bq in the vicinity of the symmetry 
axis and falls as fast as far away from the axis. This solution has wide-ranging 

applications in the literature, as, for example in the study of Kerr black holes 


22(, or 


m 


cosmology, where the possibility of the interaction of the magnetic held with the expansion 


of the universe may be considered, as for example in 


23|, |24|. 


The line element may be written in a cylindrically symmetric form (taking c=G=h=l) 


ds^ = k{r)^ dt^ — dr"^ — K{r) r‘^d(j)^ — A{r)‘^ dz‘^, (16) 


where A (r) = 1 -I- \Blr‘^ and Bq is the magnetic held. This metric rehects the curvature of 
the space-time, determined by the existence of a magnetic held in the direction. For this 
reason, a metric with axial symmetry is taken into account. 

In the limit of vanishing the magnetic held we have A(r) = 1 and the equation flT6|l 
becomes 

ds^ = dt^ — dr'^ — r'^dcj)^ — dz'^, (17) 


that is the hat Minkowski space-time written in cylindrical coordinates. So, we may dehne 


Via)(b) 


/ 1 0 0 0 \ 

0-10 0 
0 0-10 
\0 0 0 -1 / 


and choose a diagonal tetrad basis e 


(a) 


/ A (r) 0 0 0 \ 

0 A (r) 0 0 

0 0 rA(r)~^ 0 

0 0 0 A (r) y 


(18) 


( 19 ) 

















for which the equation 


is satisfied, and then it is easy to determine its inverse form 


<) = 


0 0 

A(r)"^ 0 

0 A (r) /r 


0 

0 

0 


( 20 ) 


/ A(r)"^ 0 0 0 \ 

0 
0 

\ 0 0 0 A{r)~^ J 

We are interested in studying the effect of a magnetic held Bq, that modihes the space- 
time geometry, as it is shown in the Melvin metric. If we consider a Dirac particle inside 
this held, we may also investigate the ehect of the minimal coupling in the wave equation for 
this particle considering the 4-potential In a hat space-time, a constant magnetic held 
in the z direction, Bq, that may be related to an equivalent magnetic held in the Melvin 
metric by the equation ([2]) appears if the 4-potential has the only non-vanishing component 
given by 

Observing that the term in the curved space wave equation, relative to the tetrad 

m is given by 

0 /( 1 ) 0 /( 1 ) 

( 22 ) 




equation (lT3ll becomes 
/(i) / d 
dr 


^(0) 5 ^ 


A dt A 


2A2 dr 2rA’ 


1 aA ^ 1 ^ d , y(^) d , . 1^‘^'^AA^ 


= 0 . 
(23) 


In the limit A ^ 1, this equation reduces to the usual Dirac equation for a free particle in 
a hat space-time in a cylindrical coordinate system. So, we may interpret equation (|23ll as 
a generalization of a Dirac equation for a particle inside a magnetic held, an equation that 
has been extensively studied in the literature, in many contexts, as for example 25|- 29|. 


We are interested in observing the corrections in the energy spectrum, due to the alteration 
of the geometry of the space-time, determined in eq. 

Making a transformation in equation fl23|) 

1 


= 


^rA (r) 

we obtain a simplihed form 

-y(o) Q ^(1) Q y(2)A d y^^^ d 


+ 


A dt A dr r dcj) A dz 




-I- iM -I- iq - - 


^l^ = 0, 


(24) 


(25) 


























where the usual representation for the gamma matrices is considered 



o 

o 

o 


/ 0 0 0 l\ 

^(0) = 

0 10 0 

0 0-10 

, 7^^) = 

0 0 10 

0-100 


\0 0 0 -1 y 

( 0 0 0 -A 


iy -1 0 0 0 / 
/ 0 0 1 0 \ 

^(2) = 

0 0 z 0 

0 z 0 0 

\ -z 0 0 0 / 

, 7® = 

0 0 0 -1 

-10 0 0 

\ 0 10 0 / 


(26) 


that has no dependence in the t and 0 variables. We will suppose a solution of the form 

^ = R{r) exp[—iat + ipzZ + im(()], (27) 


where a is the energy of the system, that assumes positive values for particles and negative 
values for antiparticles, pz is the momentum, m = ±1, ±2, ±3,..., a quantum number, and 
M, the electron mass. So, equation fl2^ may be written in a explicit form 

(28) 

(29) 

(30) 

(31) 


( d Tfl Q \ 

+— —I --I + (+cr — MA) Ri — PzRs — 0, 

dr r r J 

* f+W - ^3 + - MA) R2 + PzRa = 0, 

\ dr r r J 

I + ^A^A^^ R 2 + (-a - MA) R^ + pzRi = 0, 

dr r r J 

, I + !!^a2 _ ^A^A^ ) R^ + (-a - MA) i?4 - PzR 2 = 0. 
dr r r J 


In a hrst approximation we will solve the equation neglecting terms in higher orders of r, 
and then, the system of equations may be written as 

d m mB^ 2q qBq 


“’'+*1 + 7 + 2 


-r — 


Bor 


r ] i ?4 + (—(7 + M) Ri + PzRs — 0, 


d m mBl ^ 2q ^ qBo 


“ 7 “ ^ 77 ^ 7 “'" ' ^ + m)R2- PzRa = 0, 

d m mBn 2q qBo , „ , ^ ^ ^ 

—^ I +-—I-1 —r ——-—r i?2 + (~cr — M) R^ + PzRi = 0, 

dr r 2 Bor 2 

d m mB^ ^ 2q ^ qBo 


-i I +--+ —-h -—r Ri + {-a - M)Ri- PzR 2 = 0, 

dr r 2 Bor 2 


(32) 

(33) 

(34) 

(35) 






















This set of equations may be decoupled by multiplying the hrst equation by the expression 

(36) 


^ d m mBn 
dr r 2 


0 , 2g qBo 

-r + — -^- —r 


Bor ' 2 

and then, using the second and the fourth to eliminate the terms containing the spinors Ri, 
and i? 3 . The result is 

52^2 


m' {m! + 1 ) , / , 1 

^ ' ' b{ -m! + 




i ?4 (r) = 0, (37) 


dr^ ^2 ' ” ' 2 ) 4 

where b = ttiBq — qBq and m' = m — 2q/B q. In a similar way we derive equations for i?i, 
i? 2 , -R 3 that may be resumed in the form 


d^ m' {m! ± 1 ) 


-m' ± - 
2 


^ 2^2 




i?(r)=0, (38) 


dr^ V 2/4 

where the positive sign refers to R 2 (r) and -R 4 (r), and the negative one to -Ri (r) and R^ (r). 

As we may observe, equation fl38|) is similar to the Schrodinger equation, and it may be 
written as 


d'^R\ 


dr"^ 


+ {E-V,f)R{r)=Q, 


(39) 


where E = a‘^ — —pl- The term 14/ = ^ — b (—m' ± ^) + niay be identihed as 

an effective potential and as we can see, the system has the form of an isotropic harmonic 
oscillator. 

In fact, the solution of the equation fl39|l may be mapped into a 3-dimensional harmonic 
oscillator-like one in spherical coordinates. These solutions are given in terms of the associate 
Laguerre polynomials 


R{S) = 


( iViR(l/2-mV2)g-|2^1/2-m' 

/Yg^(l/2-m72)g-|/^l/2-m' 


where S = ^ and Nt are normalization constants. The energy spectrum, relative to this 
solution is 

^ = (^ +0 ^ = 0,1,2,3,..., (40) 

and using the dehnition of E given in eq. fl5^ . we obtain the energy spectrum 


a = \l + pI + 2Boq ( ^ + s ^ j -h 2B^m fn -h s -h ^ ,, 


(41) 


















where s = ±1/2 is the spin quantum number of the particle and m = 1, 2, 3..., in order to 
keep the conservation of the parity. As it was pointed before, equation fHTj) is written with 
c = G = h = 1. By making the conversion to the international system of units we have 


a = 


y/± pIc^ ± 2hc^Boq ^ 


n ± s ± ^ j ± Aneoh^GBlm f ^ ± s ± ^ ). (42) 


where Eq is the vacuum permissivity constant. In a system of units with h = c = 1 we have 


a = \l+ pI + ‘^Boq ( ± s ± ^ j ± GB^m f u ± s ± ^ , 


(43) 


were, in the last term inside the square root, the Planck mass rup appears scaling the 
magnetic held as {Bo/rripY. 


IV. RESULTS 

Studying equation fl4T|) we may recover some literature results. When the magnetic held 
and Pz goes to zero, we obtain the expression for the rest energy of the particle 


a = M(r. 


(44) 


Now, if the gravitational energy is neglected, the last term inside the square root vanishes 
and we have 

a = y/± pIc^ ± 2hc^Boq ± s ± , (45) 

that corresponds to the usual spectrum obtained from a Dirac equation in a hat space with 
a vector potential. This fact conhrms our initial considerations about eq. fl23l) . interpreting 
it as a generalization of this equation. It is also possible to analyse the situation where only 
the gravitational term is considered. 


a = \ M‘^c^ + pIc^ ± ATTEoK^GBlm ( u ± s ± 


and observing this equation it is easy to see that this term (when the magnetic held appears) 
is a negligible correction for usual magnetic helds (not so strong) for all energy levels. 

But we are interested in studying the ehect of extreme magnetic helds, as intense as the 
ones found in magnetars, of the order of B r\-/ 10^^ — [l[, |2| or the helds expected to 








be produced in heavy-ion collisions B ~ jsl-js]. The results for some energy levels, 

obtained with the exact numerical solution of eq. (|23|) (Melvin metric), considering some 
systems of interest, are shown in Table 1. Observing the Table, we can see that the effect of 
the magnetic field appears when B increases and becomes important when these fields are 
intense. 

The same calculations have been performed in a flat space-time with a minimal coupling. 
The results are shown in Table 2. As we can see, the results are essentially the same, and, 
when comparing then with the ones found in Table 1, we can only find deviations when 
considering fields of the order of i? ~ and quantum numbers as large as n ~ 10^°. 

These results are shown in Fig. 1. Observing these results, the conclusion that we may 
obtain, with a very good precision, is that the effect of the magnetic held in the metric is 
very small and may be neglected. 

As it was explained before, eq. flTTl) is an approximation for small values of r, but when 
comparing the results with the exact numeric calculations that have been performed, up 
to the precision shown in Tables 1 and 2 the conclusions that we obtain are essentially the 
same, and then we may use this equation as a very good approximation. The corrections 
from terms with higher powers of r appear beyond the precision shown in these tables. This 
fact may be seen if one observes the approximations that have been made when eq.f32)- (l35|) 
have been obtained. Terms of the order where n is an integer, appears. These 

terms increase with r^, and reaches the maximum values with some tm, that determines the 
size of the considered system. For heavy-ion collisions, for example, {BqTmY ~ 10“^^, for 
pulsars, 10“^^. The greater value obtained was for magnetars, 10“^, and the effect of these 
corrections, as it was said before, are always beyond the precision shown in the tables, so, 
we may conclude that the accuracy of the approximation is good, even for extremely large 
magnetic fields and levels. 

One remark that must be made is that when considering a covariant Dirac equation, a 
particle with gyromagnetic ratio g = 2 is taken into account. This is a good approximation 
for electrons inside stars, for example. But when considering higher energy processes, as 
high-energy collisions, deviations from this value are proposed to exist, and then, a way to 
study this question is to consider extensions of the covariant Dirac equation, as for example 
as it has been made in 2l|. This kind of approach is left for future works. 


Another aspect that may be taken into account is the Melvin metric. It is clear that 






FIG. 1. Energy levels as functions of the quantum number n considering eq. (j23p . (??) and eq. 
(j23p without a minimal coupling (solid line). 



E, (TeV) 

£iooo (TeV) 

£,(,=« (TeV) 

^Earth~0-56 G 

5.109990646783791x10“’ 

5.109991223884968x10“’ 

7691.365748802196 

Bs.„~1.5G 

5.109990650256817x10“’ 

5.109992381560156x10“’ 

13321.83625652045 

fipulsar-lO'^ G 

0.00005957925431687425 

0.001087723495543327 

3.439683331294733x10'" 

^Magnetar'lO^^ G 

0.0005957708482387581 

0.01087723376713101 

3.43968337759476x10" 


0.05957706291159353 

1.087723375512916 

3.440150976087315x10" 


Table 1: Energy levels in the Melvin metric with a minimal coupling, from the value of the 
magnetic field found in the Earth, up to the one expected to be produced in ultra-relativistic 

heavy-ion collisions. 


many of the systems that have been studied does not have the magnetic held in the form 
of the one that determines this metric, but in some regions, with intense magnetic helds, at 
least as a hrst qualitative approximation, these results must be correct, and we expect that 
more careful calculations, with the metric determined by different shapes of the magnetic 
helds, conhrm our results. 



£3 (TeV) 

£'1000 (TeV) 

£ 10 ™ (TeV) 

^Earth'“0-56 G 

5.109990646783791x10“’ 

5.109991223884968x10“’ 

7691.365748802196 

Bsu„~1.5G 

5.109990650256817x10"’ 

5.109992381560156x10"’ 

13321.83625652045 

fipulsar~IO"G 

0.00005957925431687425 

0.001087723495543327 

3.439683330827056x 10'" 

^Magnetar~10^^ G 

0.0005957708482387581 

0.01087723376713101 

3.439683330827056x10" 

Bio„s~10‘"G 

0.05957706291159353 

1.087723375512916 

3.439683330827056x10" 


Table 2: Energy levels in a flat space-time. 
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